Angular Dependence of Jitter Radiation Spectra from Small-Scale 

Magnetic Turbulence 

Sarah J. Reynolds, Sriharsha Pothapragada 
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045 

Mikhail V. Medvedev 1 
Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, 2100 

Copenhagen K, Denmark 

Received ; accepted 



1 Also: Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045 
and Institute for Nuclear Fusion, RRC "Kurchatov Institute", Moscow 123182, Russia 



- 2 - 
ABSTRACT 

Jitter radiation is produced by relativistic electrons moving in turbulent small- 
scale magnetic fields such as those produced by streaming Weibel-type insta- 
bilities at collisionless shocks in weakly magnetized media. Here we present a 
comprehensive study of the dependence of the jitter radiation spectra on the 
properties of, in general, anisotropic magnetic turbulence. We have obtained 
that the radiation spectra do reflect, to some extent, properties of the mag- 
netic field spatial distribution, yet the radiation field is anisotropic and sensitive 
to the viewing direction with respect to the field anisotropy direction. We ex- 
plore the parameter space of the magnetic field distribution and its effect on the 
radiation spectrum. Some important results include: the presence of the harder- 
than-synchrotron segment below the peak frequency at some viewing angles, the 
presence of the high-frequency power-law tail even for a monoenergetic distribu- 
tion of electrons, the dependence of the peak frequency on the field correlation 
length rather than the field strength, the strong correlation of the spectral pa- 
rameters with the viewing angle. In general, we have found that even relatively 
minor changes in the magnetic field properties can produce very significant effects 
upon the jitter radiation spectra. We consider these results to be important for 
accurate interpretation of prompt gamma-ray burst spectra and possibly other 
sources. 

Subject headings: radiation mechanisms - relativistic plasmas - gamma-rays: bursts 
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1. Introduction 

Radiation from astrophysical or laboratory plasmas depends upon the strength and 
topology of the magnetic field within the plasma, as well as on the acceleration of plasma 
particles. For relativistic plasmas, radiative output of an individual particle depends 
significantly on the particle's direction of motion, as relativistic beaming modifies the 
angular distribution of the radiation emitted to lie within a cone of opening angle I/7 
along the particle's velocity. Particles in uniform magnetic fields produce synchrotron 
radiation, which is characterized by the sweep of the radiation cone as relativistic particles 
orbit or spiral along the magnetic field lines. In turbulent magnetic fields, particles still 
radiate in a synchrotron regime as long as the scale of the magnetic field variation is longer 
than the particles' average Larmor radii. However, if the magnetic field varies on a scale 
smaller than a Larmor radius, the particles "jitter" through a series of small transverse 
accelerations without a substantial change of direction. An observer will thus be in a 
particle's radiative cone over some length of its accelerated path, until the particle direction 
is sufficiently deviated that the line of sight is outside the radiative cone. For an isotropic 
particle distribution, the rates of particles entering and exiting paths within I/7 of the line 
of sight are equal. The resulting radiation will be that of particles "jittering" with small 



randomized accelerations, and will be reflective of the magnetic fie 



sight path through the turbulent magnetic field region (jMedvedev 



d spe ctra along a line of 



2000). 



Whereas in synchrotron radiation spectra the frequencies depend upon relativistic 
beaming and the sweep of the radiation cone, in jitter radiation spectra the frequencies 
instead depend directly upon the turbulent variations of the magnetic field and relativistic 



beaming serves only to limit the o 



I/7 of the line of sight (jMedvedev 



jserv ed radiation to that emitted by trajectories within 



20001 ) . The resulting radiative spectrum differs from that 



of synchrotron radiation and, in particular, its low-frequency spectral index is not limited 
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by s = 1/3 (i.e., dW/dw oc cj 1 / 3 below the synchrotr on peak), referred to in the literature as 



the synchrotron "line of death" (jPreece. et al 



19981 ) 



Gamma-ray bursts (GRBs) are believed to be a natural astrophysical site for the 
emergence of a strong small-scale turbulent magnetic field. Relativistic shocks generate 



strong magnetic fields that are r andom on very sma 



(particle streaming) instability (jMedvedev fc Loeb 



(sub-Larmor) scales via the Wiebel-like 



19991 ) . In such an instability, a small 



magnetic field perturbation in the plane transverse to the motions of counter-streaming 
particles results in the development and growth of filamentary current structures up to 
a saturation point at which they may persist for some time longer than the dynamical 
time-scale of the system, creating an extended region of small-scale tur bulent magn e tic fie ld. 



This has been st u died extensively in nume r ical P I C simulations (e.g . , 



Nishi 



cawa. et al 



(2005); 



(2003); 



Chang et al. 



Frederiksen et 



(2008) 



al 



J200J); 



Medvedev. et al. 



Silva, 



(120051 ): 



et al. 



(2003); 



Spitkovsky 



Spitkovskyi (120071 )) of both baryonic and pair plasmas. (See 



Medvedev fc Spitkovsky! (120091 ) for a discussion of the applicability and implications of such 



simulations for GRB physics.) Recent simulations of magnetic reconnection in pair plasmas 
have also shown the generation of these strong small-scale magnetic fields via the Weibel 



i nstability acting on t 



(jSwisdak. et al. 



le streams of accelerated particles in the reconnection exhaust funnels 



(I2008h : 



Zenitani &; Hessd (120081 ) ) . We speculate that magnetic reconnection 



events in a GRB may produce electron-positron plasmas in situ, even in initially lepton-poor 
plasmas. For the resulting regions of strong small-scale magnetic field, synchrotron theory 
is inapplicable and jitter radiation theory must be considered. Thus, regardless of the 
particular model of a GRB (e.g., baryonic, leptonic, magnetic, etc.), the jitter radiation 
mechanism coupled with the relativistic kinematics of the ejected material represents a 
viable phenomenological model. It has recently been shown that s uch a model reproduc es a 



number of spectral features of GRB light-curves remarkably well (jMedvedev et al 



20091 ). 
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Following the approach of 



Medvedevl (120061 ). we define a geometry where the local 



filamentation axis (the local axis along which counter-streaming particle motion occurs) 
lies in the z-direction so that the field perturbation from the shock or reconnection event 
is amplified by the Weibel instability in the xy-p\a.ne (this might be at or upstream from 
a shock front lying in the xy-plane and propagating in z, or upstream from reconnection 
exhaust funnels with filaments pointing in z direction and the Weibel fields being 
perpendicular to it, as suggested by PIC simulations). The resulting amplified magnetic 
field is randomly oriented i n the x?/-plane and i ndepe ndently generated at each position 



in z (sh own theoretically in 



such as 



Nishikawa. et al. 



(120031 ) . 



Silva. et al 



Medyedev fc Loebl (119991) an d confirmed in P 



(120031 ). and 



Frederiksen et al. 



C simulations 



(12004m . The 



decoupled behavior of the magnetic fields along the filamentation axis (z) and in the plane 
(xy) transverse to it means that the resulting Fourier spectra are independent of one 
another and the overall field distribution is highly anisotropic. Qualitatively, the spectrum 
of the Weibel-generated magnetic field in the direction transverse to the filamentation 
axis has been shown to rise and then drop at a scale of order the plasma skin depth 



( IFrederiksen et al 



20041 ). The magnetic field thus has a general spectral form that may be 



parameterized as 



k 



2a 



where k± = (k 2 . + k 2 ) 1 / 2 , and a± > and (3± > 0. In this form k± refers to the magnetic 
field wavenumber in the plane transverse to the filamentation axis and k±, a±, and /3± 
are free parameters controlling the spectral break and the soft and hard spectral indices, 
respectively. The field along the filamentation axis is in general unknown but we expect it 
to be of a similar form, with independent free parameters k\\, a\\, and f3\\: 



2a II 



(2) 



where again a\\ > and (3\\ > 0. A plot of f z (k\\) is shown for a particular choice of 
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parameters in Figure 1(a) and demonstrates the basic behavior of this function. 



The variables k and k are presumed to be unitless, their units fco havi ng been separated 



o a 

into a normalizable coefficient k . We have modified this form from the 



Medvedevi (12009 ) 



paper, in which the power in the denominator was simply (3 and it was required that (3 > a. 
In these spectral forms used here, the asymptotes of the functions given in equations [T]r2] 
are: 

/(*) = < (3) 



if k « k, 
if k >> k. 



by 



The theo r y of U tter radiation with the above magnetic field spectra has been utilized 



Medvedevi (120061 ) to derive the basic equations describing the jitter spectrum and 



demonstrate its dependence on the angle 9 between the viewing angle and the local 
filamentation axis. Taking i nto account our m odification of the form of the magnetic field 



spectra, the analytical work (IMedvedev 



20061 ) indicates that the jitter radition F u spectrum 



should have the following general properties: 

1. two breaks, with locations depending on k±, k», and 9, 

2. a high-energy spectral index /3', where f3' approaches (3\\ as 9 goes to and (3' 
approaches as 9 goes to tt/2, 

3. a low-energy spectral index a', where a' approaches 1 as 9 goes to and a' approaches 
as 9 goes to tt/2. 



The resulting radiation spectra from numerical cal culations fo r a lim ited selection of 



Medvedevi (120061 ). Here we present 



parameter choices and viewing angles were presented in 
the results of more extensive numerical calculations illustrating the full range of jitter 
spectral variation due to viewing angle for a chosen set of field parameters, followed by a 



more thorough exploration of the magnetic field parameter space and its implications for 
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the jitter radiation spectrum. Section II of our paper presents the calculated acceleration 
spectra for a range of viewing angles 9. In Section III we develop the connection between 
the acceleration and radiation spectra through an analytical treatment of a simple linearized 
model of the acceleration spectrum. In Section IV we present the radiation spectra 
calculated for various viewing angles 9 and analyze the spectral variation with viewing angle 
by means of a five-parameter spectral fit. In Section VI we show the effect of variations in 
the magnetic field parameters upon the radiation spectra. Section VII presents discussion 
and conclusions. 



2. Acceleration Spectra 



The equations for calculating the spectrum of a particle's acceleration due to magnetic 
field turbu l ence generated by a relativistic Weibel-type instability were developed in 



Medvedevl ( 120061 ). but only calculated in full for a single representative oblique viewing 
angle of 9 = 10° in between the head-on (9 = 0) and edge-on (9 = 90°) cases. Here we 
calculate the acceleration spectra for a more complete range of intermediate viewing angles 
in order to explore the spectral progression with 9. 



As derived in iMedvedevI ( 120061 ). the volume-averaged temporal Fourier component of 
a particle's acceleration due to the Lorentz force for the static magnetic field generated by 
the relativistic Weibel instability is 



(|w^| 2 ) = {2<kV)~ x J \w k \ 2 5(uf + k- v)dk 



C 
2^ 



(1 + cos 2 9) J f z (k ]{ )f X y(k ± )6{uj/ + \i-v) d 2 k L dk\\ 



(4) 
(5) 



iui't 



where k is the magnetic wavevector, v is the particle's velocity vector, and = J we 
For the case of a shock viewed at an oblique angle 9 from the normal of the shock plane 
(defined as the 2-axis), k • v = k x v sin 9 + k z v cos 9, where we have defined the x-axis so that 
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the velocity vector v lies in the xz-plane. We can then use the delta function to substitute 
for either k x or k z . This becomes: 

W -'| 2 > = o i ° .i ( 1 + cos2g ) / fz (-^-a+k^tane) f xy ((k 2 x + k^ 2 )dk x dk y (6) 
x ' 2ir\v costal J \vcos9 J u 



or 



C f / / / Uj' \ 2 \ 1 / 2 \ 

(|w^| 2 ) = - | . ,. (l + cos 2 fl) / f z (k z )f xy (((—- + k g cote) + k 2 y ) )dk y dk z (7) 

where C is an arbitrary normalization constant, proportional to {B 2 ). The two forms are 
equivalent; however, as we approach 9 = or 9 = tt /2, the calculation is more convenient if 
one avoids denominators approaching zero by choosing the appropriate form. It should be 



noted that neither form is valid 
separately as in 



Medvedevl ((20061). 



or the endpoints 9 = or 9 = tt/2, which must be treated 



For a single radiating particle, we plug equations (jTJ) and (j2J) for f xy and f z into 
equations and (JTj) to obtain: 

where 



2t« 7-c i-oo («| cot 2 + (^5 + fc r )2)"n+< 3 . («» + *i + tj)^ 



COS0 



or alternatively 



(10) 



where 



= (oot g )-^(l + co^) 
sin 



This may then be numerically integrated to produce the acceleration spectrum as a function 
of u'. 
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We normalize the wave- vector to a dimensional constant fco, hence k is dimensionless, 
as discussed above. The spectral form remains the same with k = n'ko and a normalizable 
factor A)q in front. The frequencies u are normalized to uj such that v = u>o/k = 1. 
The parameters oti, and /c» in our magnetic field spectra are in general unknown free 
parameters that depend upon the spatial distribution of the magnetic field in the turbulent 
region, and the parameters Kj, corresponding to the dimension correlation lengths of the 
magnetic field in the turbulent region, will vary even within an astrophysical system and 
affect the resulting peak energy of the radiation spectrum. For our initial calculations of 
spectral variations with 9 we choose reasonable values of the parameters as follows: 

K± = W || = 10, (12) 
a± = a\\ = 2, (13) 
P± = fl| = 1-5; (14) 

These parameters are varied (jointly and individually) in Section [5] to explore the resulting 
effect on the calculated radiation spectrum. 



We present the results for varying 9 in Figure 2(a) Figure 2(b) shows further detail 
of the critical region around the peak of the acceleration spectra. The graphs show 
linearly-connected log e — log e data, with data point intervals of 0.1. We have arbitrarily 
normalized the spectra so that the low- A; part of the 9 = 2 spectrum asymptotes at unity. 
The acceleration spectra are flat for low a/, then at a certain u' they turn rapidly into a 
sloped region, then there is a second transition to a steep decline for high uj' . For 9 near 
the spectra have a clear peak, but as 9 increases the peak recedes and eventually disappears 
altogether. As 9 approaches tt/2 a peak again becomes evident but the location of the peak 
has shifted by about 0.4 from its position for low 9. The transition point for flattening at 
low uj' appears to move off rapidly to lower to' as 9 approaches 0. 

Figure [3] plots the amplitude of the acceleration spectra at our lower calculation 
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boundary, at the approximate location of the peak for small 6, and at the approximate 
location of the peak for 9 close to it/2. The crossing of the lines on this graph correspond to 
spectral transitions as the peak disappears and then eventually reappears at a new location 
for higher 9. The graph of the slope of the log — log plot in Figure H] also illustrates the 
disappearance and reemergence of the peak, but further shows that even for the unpeaked 
spectra there is a flattening of the slope of the mid-range 9 spectra in the region in between 
the positions of the peaks that appear at higher and lower 9. For unpeaked spectra, there 
remains a transition region of some extent between the low-frequency and high-frequency 
power laws; consequently, the unpeaked spectra may still be better fit by division into three 
power law regions as opposed to two. 

An analysis of equations (jHJ) and ( flOl) indicates that u' functions as a shift of the 
center of the magnetic spectral form in which we have made the delta-function substitution. 
Since this uj' factor is always positive, the function's offset always occurs in the direction of 
negative wavenumber components k x or k z . An example of this behavior for the product of 
two functions such as equations (0Q) and (J2J) over a range of offsets is shown in Figure EJ 
The resulting integral is highly sensitive to the shape of the two functions and the offset 
between them, which control the resonance like behavior that produces the resulting peak 
and transition points in the acceleration spectrum. The angle 9 plays a role in determining 
both the width and the offset of one function relative to the other: the k| cot 2 9 or hi\_ tan 2 9 
terms in the denominator influence the width of the function under consideration, while the 
offset is given by - D or w . a . These are linked such that as 9 increases, both the width 

° J v cos v sin t) ' 

and the offset of the function containing 9 increases. The spectral indices a±, a\\, f3±, and 
/3n can also influence the width of the functions and hence also affect the location of the 
transition points in the acceleration spectra. The overall effect of such variations on the 
resulting jitter radiation spectrum will be explored more in Section 5. 
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Summarizing the results of this section, we find that the acceleration spectra are 
generally characterizable to a good approximation by division into three regions: a flat 
low-/c region and two power law regions, as shown in Figure These may altogether be 
described via one amplitude, two non-zero spectral indices being functions of as and /3s of 
the field spectra, and two breaks, which depend primarily on the magnetic spectral peaks 
k±, K\\, and the viewing angle 9. In the next section we develop a simple linearized model 
of the acceleration spectrum using these five parameters and show that they can be used to 
analytically predict the behavior of the final jitter radiation spectrum. 



3. From Acceleration to Radiation Spectra 



The angle-averaged emissivity of a relativistic particle undergoing a series of small 
transverse accelerations not substantially affecting its overall velocity is as follows 



( ILandau &: Lifshitz 



1971 



Medvedev 



dW 
duj 



e 2 oj 



2000L and others): 

2 



W,„/ 



27TC 2 7u,/2 7 2 ^ 2 



1 - 



UJ 



L) 



a/7 2 



2a/ 2 7 4 



du' 



(15) 



Analysis of the high-u; and low-u; asymptotic behavior o 
viewed head-on and edge-on, as carried out in 



Medvedev 



the r adiation spectra for shocks 



(120061 ) . yields: 



for 9 = 0, 



for 9 = n/2, 



dW 

duj 



oc 



UJ 1 

uj~ 2 ^ 



if to <C /t||t>7 2 
if uj ^> /t||t>7 2 



dW 

duj 



,0 



uj" if u <C /tj_f7 2 

u -2P±+i if a; > k ± vj 2 



(16) 



wher e 9 is the angle between the line of sight and the normal to the shock front (jMedvedev 
20061 ) and a±, ac\\ > 1/2. We thus expect the radiation spectra at oblique viewing angles 
to vary between these two forms, dominated by the parallel or the transverse spectra as we 
vary between the two extremes. 
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The u dependence contributed by the integral primarily originates in the lower limit 
lj/(2'j 2 ) and where it falls on the acceleration spectrum. In section [2], we found that the 
acceleration spectrum could be simply characterized as three regions of approximately 
power-law behavior: a flat initial amplitude at low logo;' (Region I), a region of positive or 
negative slope (Region II), and a region with a more steeply negative slope (Region III). 
Using this to make a simple approximation for our acceleration spectrum in three regions, 
we can calculate an approximate analytical solution to equation [15j 

We define our simple approximation to the acceleration spectrum using five free 
parameters, all of which may be fit to the spectrum on a log e — log e plot: j? = loga is 
the amplitude of the low-u/ limit; 1\ = logu^ is the transition point between the first and 
second regions; Si is the spectral index in the second region; T2 = logtu^ is the transition 
point between the second and third regions; and —S2 is the spectral index in the third 
region (52 > 0). The acceleration spectrum then has the form (shown in Figure [6]): 

for J < e Tl : 

<|w^| 2 > = e* =a , (17) 
for e Tl < u' < e' T2 : 

(M 2 )=e*-**u* = aa(^Y , (18) 



for uj' > e T2 : 



Ow^i 2 ) = e ^™+, 2 .r V -, 2 = ao My 1 f^_y S2 (19) 

In Figure [7] we show example models for the acceleration spectra at 9 = 10 degrees and 9 = 



60 degrees and their comparison to the original acceleration spectra (as in Figure 2(a)) for 
a particular choice of fitting rules. 

To first order (neglecting the second and third terms in Eq. (fl5|) ). the resulting 
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radiation spectrum is as follows: 



for cj/27 2 < e Tl 
'dW 



for e' Tl < cj/27 2 < e r2 
'dW 



du 



11 



27TC 2 

ape 2 
2tvc 2 



2ttc 2 

age 2 
2vrc 2 



(2 7 2 ) - u 



2 7 2 - 51 



+ 



(5i + 5 2 )e 



-Tl5l+T2(5l-1) 



(5i-l)(5 2 + l) 



5l +5 2 



Si - 1 V^i/ (5l - l)(s 2 + i) Wl 



LO 



I \ Si 



LO 



7 m 



-lo 



Si 



.(5i-l)(2 7 2 )^- 1 
(2 7 2 )^ 



LO 



(51-1) w 



,5i - 1 5 2 + 1 
5i +5 2 



,0-2(^1-1) 



(5i-l)(5 2 +l) W 



CUo \ to 



LOU 



, (21) 



for cj/2 7 2 > e T2 : 



dto 



A-?lSi +1l{Sl +52) 



(27 



2U+52 



J/J 



27TC 2 

a e 2 (2 7 2 ) 1+ ^ 
2ttc 2 |5 2 + 1| 



5 2 + l 



-a; 



-J2 



to 



(22) 



The radiative power spectrum (equations ( 1201) - (1221) ) obtained analytically by our 
fit-based approximation for the acceleration spectrum agrees with that obtained via full 
numerical integration in the following section within about 10%. We have chosen our 
fitting method for the acceleration spectrum to most closely capture the spectral indices; a 
different choice of fit may allow for a better determination of peak positions. 

Calculations of such spectra for parameters fitted to our acceleration spectra at 9 = 10 



degrees and 6 = 60 degrees are shown in Figures 8(a) and 8(b) The region boundaries in 



these figures indicate that the radiation spectrum cannot be described by a simple linear 
approximation in the three regions that were defined by the breaks in our acceleration 
spectrum. The key features (position of spectral peak and spectral breaks) of the radiation 
spectrum originate in the additional terms in Equations (|20|) and (I2TI) . Consequently, the 
transition points in our radiation spectrum do not directly correspond to the transition 
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points in the acceleration spectrum. 

The asymptotic behavior of the spectrum at high and low energies can be easily 
obtained from Eqs. ([20]) and ([22]): 



for cj/27 2 < uj[ : 



1 

2 



duo 

for cj/27 2 > oj' 2 : 
' dW 
duj 



oc u°, (23) 



oc oj- S2 (24) 



in 

Thus the high and low-energy asymptotic behavior of the radiation spectrum will be 
identical to the high and low-energy behavior of the acceleration spectrum. 

The behavior of the spectrum in the intermediary Region II can be solved for as well. 
We take the derivative of Equation (]2~TT) to solve for the position P = \og(u p /2'j 2 ) of the 
spectral peak, in cases where it exists: 

"5i(5 2 + l) 



Si 



log 



S1 + S2 



+ ?2 (25) 



We note that the peak position becomes undefined in the case 5i < 0, \si\ > S 2 , for which 
Equation (T2~T]) is everywhere decreasing. The case s± < 0, \si\ < 5 2 , in which the acceleration 
spectrum would decline more steeply in Region II than in Region III is impermissible, so we 
find that the radiation spectrum will be unpeaked whenever the mid-range spectral index 
Si of the acceleration spectrum is negative (i.e., for 5i < 0). 

An exploration of the behavior of peak point P relative to the region boundaries shows 
that for certain values of 9 the peak of the Region II function exists, but has crossed 
the boundary into Region I and consequently does not appear in the resulting radiation 
spectrum. Thus, while a calculation of P from the results of fitting the acceleration 
spectrum appears to indicate the re-emergence of a peak in the radiation spectrum as 9 
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approaches tt/2, this peak falls beyond the Region II lower boundary and is not observed. 
We note that the first term in Equation (1251) is negative for both < Si < 1 and s% > 1, so 
the peak is always located below the transition point between Regions II and III. 

An analysis of the behavior of Equation (12TT) below the Region II peak indicates the 
following behavior: 

dW\ f^ 1 if ^i > 1, 



for lo[ < lo/2Y < uj' 2 : — oc { (26) 

\ duJ J ii [lo Si if5i<l, 

Thus, from Equations f[2"5j) - ([2"6"j) . it is evident that the spectral indices 5i and 52 of 
the acceleration spectra will generally correspond to spectral indices Si and S2 in two 
power law regions of the radiation spectra. In the case of the high-frequency spectral 
index S2 this correspondence is exact; however, the relation between mid-range (i.e. 
intermediate-frequency) spectral indices Si and s± is modified an upper limit of unity on si 
and also breaks down when the first term in equation [25] is undefined or larger in magnitude 
than the distance between the acceleration's spectral transition points 1% — T\. 

The asymptotic form in Equation ( l23i) suggests that we may neglect the second and 
third terms in Equation (120]) and solve for the transition point T = log(o;i/27 2 ) at which 
the dominating term in Region II becomes significant. We find that: 

for Si < 1 : 

=(i--5iK 1 



T =!]og(l (27) 
Si 

for Si > 1 : 

(-51-1)0*2 + 1) { < 



2 



(5i-l)(5 2 + l 



log 



+ 5iT 1 + (5i-l)T 2 (28) 



Si +5 2 

We have demonstrated that a simple approximation for our acceleration spectrum 
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allows us to analytically obtain some of the key features of the resulting jitter radiation 
spectrum, notably that it will have a similar three-region form with dW/ du oc uj° for small 
uj and dW/du oc oj~ S2 for large u and a possibly-peaked transition region. Unlike the 
acceleration spectrum, the intermediary region in the jitter radiation spectrum will have a 
maximum slope of 1 and may be unpeaked at angles 9 at which the acceleration spectrum 
was peaked. We emphasize that our spectral calculations are all for a single emitting 
electron, not a power-law distribution of electrons. Yet, the power- law photon spectrum 
emerges at high energies (above the second jitter spectral break oo 2 ) in this jitter mechanism, 
in contrast to the synchrotron exponential spectral decay above the synchrotron frequency. 



4. Analysis of Radiation Spectra 

Now we turn to full numerical calculations of the jitter radiation spectrum, generated 
by successive numerical integrations of equations (IT5|) and (JSJ) or ffTUj) The results for 



varying 9 are presented in Figure 9(a), with data point intervals of 0.2 on the logo; scale. 



Once again we have normalized the spectra such that the low-energy asymptotic value of 



the 9 = 2 spectrum is unity. The detailed view of the peak region in Figure 9(b) shows 
results in this region for intervals of 0.05 in log(o;). The spectral shapes and trends are, as 
expected, much like the acceleration spectra but broadened and flattened overall. No peak 
reemerges in the spectra as 9 approaches 7r/2. Figure [TUl shows the vF v spectrum such as is 
commonly presented for GRBs and used in GRB spectral analysis. 

These jitter radiation results show a significant evolution in the spectrum emitted 
at different viewing angles relative to the main filamentation axis of the magnetic field 
spectrum. Note that this viewing angle effect in our calculations is entirely due to particles 
with velocities directed along the line of sight providing the dominant contribution to 
the radiation emitted to any particular viewing angle. We are neglecting the angular 
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distribution of the radiation emitted by each particle and using the angle-averaged 
emissivity for the spectrum emitted in the forward direction by particles with a particular 
orientation angle relative to the magnetic field filamentation. 

Like the acceleration spectrum, the radiation spectrum can be generally described 
in terms of three regions (two spectral breaks), and an amplitude or slope in each. To 
conveniently summarize the spectral features and their evolution, we have developed a 
five-parameter fit which describes the spectral behavior in these three regions, which we 
designate as Regions R-I, R-II, and R-III to avoid confusion with the Regions I, II, and III 
as defined earlier for the acceleration spectrum. (Recall that the transition points for the 
acceleration spectra do not correspond to the apparent transition points in the radiation 
spectra.) As before, we have chosen our technique to optimize our results for the spectral 
indices, rather than the spectral transition points. 

Since for jitter radiation the soft spectral index varies continuously and approaches 
for low ou, the results of a simple two-region fit to these spectra would depend significantly 
upon where the lower bound of the data window falls relative to the peak. Within a 
fixed data window, a two-region spectral fit would tend to produce an artificial reduction 
in the soft spectral index for spectra with higher-frequency spectral peaks or breaks. 
Unfortunately there is no simple way to characterize the behavior of the middle range of 
the spectrum because of the transition from peaked to unpeaked spectra as 6 varies. Even 
in unpeaked spectra, the extent and curve of the transition region between the flat low 
uj' part of the spectrum and the strongly negatively-sloped high uj' part of the spectrum 
varies substantially. Consequently, we have chosen to model the unpeaked spectra still as a 
three-region spectra rather than solely by its upper and lower asymptotes. We characterize 
the spectra by defining three lines (each requiring a slope and a reference point) and finding 
the transition points at which they intersect. 
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Region R-I (flat, amplitude A): The low-frequency region R-I is flat, with a slope 
close to zero. To describe this region, we take our initial calculated amplitude to be A, 
the low-frequency amplitude. For very small 9 our lower calculation boundary may be 
insufficient to capture the initial flat part of the spectrum, since our first spectral transition 
point approaches -oo as 6 goes to 0. 

Region R-II (positive or negative slope si): The intermediate-frequency region R-II may 
have positive slope resulting in a peak or a slight negative slope (of notably less magnitude 
than the slope in region R-III). Since not all the spectra are peaked, we have chosen to 
avoid using the peak value for our fit. Instead we define the "drop point" as the region 
where the second derivative reaches its minimum value. As the place of largest negative 
change in slope, this coincides well with the "knee" or second break of the function, and 
is always at slightly higher frequencies than the peak itself. We then find the slope and a 
reference point in this region by either: 

• Method a: for peaked spectra, we take the maximum value of the numerical derivative 
and its associated data point. 

• Method b: for unpeaked spectra, we take the average value of the numerical derivative 
in the region between the drop point and the "deviation point" where the spectrum 
first drops below A — 0.01 (this corresponds to a deviation of about 1% from its 
original value). Our reference point is the data point halfway between or next highest 
to halfway between the deviation and drop points. 

Region R-III (negative slope, defined as -s 2 ): The high-energy frequency region R-III 
has a large negative slope compared to the rest of the function. This slope is still changing 
over the region close to the second spectral break that we are considering, so we determine 
a representative slope by calculating the slope of a line between the drop point, which is 
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the minimum of the numerical second derivative, and the higher frequency point that is 
the absolute minimum of the numerical second derivative (the data point closest to where 
the second derivative crosses zero in this region). These points are well-defined for all our 
radiation spectra as long as the calculation boundary extends a couple orders of magnitude 
in e above the drop point. Either point may be used as a reference point in this region. 

The first spectral transition point T\ is obtained by solving for the intersection of the 
lines defined in Regions R-I and R-II, and the second spectral transition point T2 is obtained 
by solving for the intersection of the lines in Regions R-II and R-III. 

We have chosen to work with the F v spectrum because of the convenience of its 
distinctive flat (spectral index of 0) initial amplitude at very low frequencies, but it is easy 
enough to translate F v spectral features into spectral features of the vF v spectrum or the 
photon spectrum N(E), as the spectral indices will simply be increased or decreased by 
1 and the transition points between the power law regions will roughly coincide, with a 
slight shift ba sed on normalizat ion. In terms of the Band function fit commonly used for 
GRB spectra ( Band et al. 199^ \. the relation between the high-energy spectral indices is 



f^Band = S2 — 1. The relation between the low-energy spectral indices is complicated by the 
fact that the Band function is a two-region fit and not sensitive to multiple spectral indices 
below the spectral peak; consequently asand will range between -1 and s\ — 1 depending 
on where the data fitting window falls relative to our first spectral break t\. The vF v 
peak energy, which is E p in the Band function, will correspond to a slope of -1 in the F v 
spectrum, and will lie roughly in the vicinity of the second spectral break T2 in our fit. 

Figures [12] - [14] show spectral fit results obtained using our above technique on the 
radiation spectra. We have also applied the same technique to the acceleration spectra 
presented in section [2] and plotted them for comparison. In addition to data resolution 
effects, some discontinuity in fitting the peaked vs. unpeaked form of the spectrum is 
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unavoidable and is reflected in our results. 

Figures d2] and indicate that the amplitude A and the high-frequency spectral 
index S2 are close in both their values and their evolution with 9 for the two types of 
spectra. The mid-range spectral index s\ varies similarly in Fig. [13] for both spectra, but 
appears to approach different asymptotic values as it approaches 9 = and 9 = n/2, clearly 
showing the expected S\ < 1 limiting behavior. The second spectral break T2 shows that the 
radiation transition point tends to be about 1/2 a power of e lower than the second spectral 
behavior for the acceleration case, but shows similar evolution with 9 in both cases. 

Figure [T7] shows the angular dependence of the spectral peak in both our acceleration 
and radiation F u spectra, the vF v spectral peak (peak data point in Figure [10] and the 
drop point, which we have defined as the minimum in the numerical second derivative. 
We clearly see the usefulness of the drop point in tracking the spectral behavior across 
the full range of 9, and that it closely tracks the behavior of the vF v spectral peak (E p ). 
Both Figures [T7] and [18] clearly show the re-emergence of peaked acceleration spectra as 9 
approaches 7r/2 and the lack of peaked radiation spectra for similar values of 9. 



5. Spectral Features and Exploration of the Spectral Parameter Space 

We have explored the influence of changes in the magnetic field spectral parameters 
on the acceleration experienced by the particle and hence its resulting radiative profile. 
Sections [2] and [4] presented the acceleration and radiation spectra calculated from magnetic 
field spectra of the form given in Equations ([1]) and (J2J) with our original choice of 
parameters a = a± = a\\ = 2.0, j3 = (3± = (3\\ = 1.5, k = k± = k\\ = 10. In this section, we 
present the results of varying these parameters. We vary the joint parameters a, /3, and k 
in the parallel and perpendicular magnetic field spectra. We also vary the parameters an, 
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a±, and individually. Finally, we vary the ratio X = k±/k«. For each variation of 
the initial parameters, we have calculated the radiation spectrum for three representative 
angles at 9 = 10°, 45°, and 80°. The results are presented according to their impact on 
the characteristics of the radiation spectrum as developed in section m namely the initial 
amplitude A, the spectral breaks n and r 2 , and the spectral indices si and $2- We also 
present the results for the peak "strength", the height of the spectral peak above the initial 
amplitude. 

The spectra are divided by (B 2 ) oc J f(n\\)f(K_ L )d 3 k oc k 2 ^ ||+/3± - ) to appropriately 
normalize the amplitudes relative to one another, but in all cases we have arbitrarily 
normalized the final spectra such that the low-energy asymptotic value of the 9 = 10 
spectrum with our original choice of parameters is unity (zero on the logarithmic scale). 

As in Section HI the initial, low-frequency amplitude A is the first calculated value of 
the angle-averaged radiative power emitted per frequency dW/du. This value is generally 
a good approximation for the asymptotic value of the function as it approaches lower 
uj, though it may deviate somewhat from this value for 9 approaching 0, as the spectra 
becomes sloped rather than flattened at our lower calculation boundary in uj. 



Among the resulting figures 19(a) |2"UI variations in the magnetic field parameters 
produce the largest effect on the low-frequency amplitude, causing changes of about 4 
orders of magnitude in e when varied individually via the ratio K = k±/k», and up to 7 
orders of magnitude in e when varied together Variation with changes 

in the magnetic field spectral indices a.i and fa are small in comparison, on the scale of 
about 1-2 orders of magnitude. The amplitude increases with increasing K for 9 = 10° and 
generally decreases with increasing K for 9 = 80°; thus, it increases when k± dominates at 
small 9 and when k\\ dominates at large 9. 

The mid-range spectral index in peaked spectra is the maximum slope below the 
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peak, as determined by taking the numerical first derivative of our calculated values. For 
unpeaked spectra, we find the spectral index as the average slope between the point at 
which the graph deviates by more than 0.01 from the initial amplitude A and the drop 
point at which the numerical second derivative reaches a minimum. 



Figures 21(a) 21(c) present the effect of magnetic field parameter variations on s\. As 
can be seen in Figure [22], the mid-range spectral index is strongly affected by variations in 
both the parameters especially for 9 = 10°, and in the relative strength of the Ki values. 
The peaked 9 = 10° spectra is notably more sensitive to magnetic field variations than 
the unpeaked 45° and 80° spectra. The ratio K = k±/ k» has the largest influence at all 
three representative viewing angles, with si increasing as Ku > We note that even at 9 
approaching it/ 2, we obtain a positive slope (and hence a peaked spectrum) for K = 1/10. 

For both peaked and unpeaked spectra, the high-frequency spectral index -s 2 is 
determined by taking the slope between the drop point and the absolute minimum of the 
second derivative above the peak (i.e., the closest data point to where the second derivative 



crosses zero). Figures 23(a) - 23(c) show the effects of variations in the magnetic field 
parameters on -s 2 . We find that as expected analytically, this spectral index is primarily 
influenced by the magnetic field parameters In particular, s 2 is most strongly influenced 
by /3±, the high-wavenumber spectral index of the magnetic field spectrum transverse to the 



current filamentation. As seen in figure 23(b), (3\\ affects s 2 only at small angles 9, and its 



influence even then is less than that of varying (3±. The apparently strong influence of k is 
largely an artificial effect as the k parameter's strong shifting of the function (as indicated 
in our analysis of the spectral breaks below) towards higher frequencies interferes with the 
calculation of s 2 by shifting the absolute minumum of the second derivative outside our 
calculation boundaries. This causes an artificial reduction in the steepness of the slope for 



k = 100, as evident also in Figure 23(c 
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We calculate the first spectral break (i.e. transition point) as the intersection between 
the line log(|w^/| 2 ) = A and the line of slope s\ through the point of maximum positive 
slope for peaked spectra or through the data point in the middle of the range over which we 
averaged to obtain slope s± for unpeaked spectra. (If the middle of the range does not fall 



on a data point, we take the next larger data point.) We find, as shown in figures 25(a) 126} 
that the first spectral break is strongly influenced by changes in a and k or the K-ratio 
K. The break position shifts to higher frequency by about an order of magnitude in e as 
we increase a from 1 to 10 jointly in the parallel and perpendicular magnetic field spectra. 
In varying the oti separately we see that ot\\ has a larger influence at 9 = 10° and a± has a 
larger influence at 9 = 80°. Increasing the k jointly by powers of 10 results in shifting the 
first spectral break to higher frequencies by roughly 4 orders of magnitude in e. Varying 
the k parameters relative to one another results in a similarly strong shift, towards higher 
frequencies for k± > k\\. 

The second spectral break r 2 is calculated as the intersection between a line of slope s± 
through the point of maximum positive slope (for peaked spectra) or through the mid-point 
of the averaging region (for unpeaked spectra), and the line of slope — s 2 through the "drop 
point" at which the second derivative reaches a minimum (i.e., the largest negative change 



in the slope). Our results (in figures 27(a) 28(c)) indicate that the second transition point 



is most strongly influenced by the Ki varied jointly or via the ratio K. The low-wavenumber 
magnetic field spectral index also demonstrates a fairly strong influence, with larger a 
shifting r 2 to higher frequencies. A comparison of the influence of a on the two spectral 



break points (as seen in Figures 25(a) and 27(a) ) indicates a very similar shift in both break 



points; thus increasing a shifts the entire spectrum towards higher frequencies. 



Figures 29(a) 29(c) show the variation in the peak strength (which we have defined as 



the height of the spectral peak above the low-frequency amplitude A with changes in the 
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magnetic field parameters. For = 10°, the only of our three representative angles that has 
peaked spectrum for K — 1, peak strength increases with increasing a and j3. Individually, 
increasing fl± has the largest effect in increasing the peak strength, while increasing a± 
lowers it. Similarly, increasing an increases the peak strength while increasing /3ii lowers it. 
The largest effect overall is produced by variation of the ratio K between the perpendicular 
and parallel field parameters For k\\ > k±_ (i.e. K < 1), the peak strength appears 
to persist to higher angles 0, while for k,j_ > ku the peak can be small or non-existent 
even at = 10°. Thus the ratio between k± and k», the respective peaks of the magnetic 
field perpendicular and parallel spectra, strongly influences the progression of the spectral 
evolution between its = and = 7r/2 limiting values, as expected from our earlier 
analysis in Section [2j 

We have seen that relatively minor changes in the magnetic field spectra can produce 
very significant effects upon the jitter radiation spectra, particularly in the appearance of 
the spectral peak or break region. Furthermore, while we have included in this section only 
spectra from a few representative viewing angles 0, the angular dependence demonstrated 
suggests that the connection of such features to the transverse or parallel magnetic field 
spectra can be tested by observing their variation with viewing angle. 

6. Conclusions 

We have calculated the angle-averaged power spectra of jitter radiation emitted by 
a single relativistic electron undergoing small Lorentz-force accelerations transverse to its 
overall velocity. Note that the obtained spectra are equivalent to the ensemble-averaged 
spectra per one electron from a collection of monoenergetic relativistic electrons. The 
radiation spectra are calculated using a smoothly connected broken power-law model 
of a magnetic field mimicking the structure of magnetic fields generated by the Weibel 
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instability. The shapes of the resulting jitter radiation spectra are shown to depend on the 
magnetic field spatial spectrum and to vary with the angle 9 of the electron velocity (being 
also the line of sight) with respect to the direction of the field anisotropy (z-axis). The 
effect of varying parameters in the magnetic field spectra has been explored and indicates 
that the jitter radiation spectral features, such as the strength of the spectral peak or the 
extent of a sloped transition region, are quite sensitive to the parameters controlling the 
magnetic field spectra. 

Despite the high sensitivity of the jitter radiation spectra to the magnetic field sp atial 



2006|), one 



spectrum or, in general, the field correlation tensor, i^-(k) = B^B^ (jMedvedevI 
can draw some fairly robust conclusions. When the parallel and perpendicular magnetic 
field spectra are similar, one has just four essential parameters, their low-/c and high-fc 
spectral slopes, a > 1/2 and j3 > 0, the peak representing a typical correlation length k, 
and the viewing angle of the line of sight with respect to the magnetic filament direction, 
9. The power (i.e., F v ) spectrum produced by monoenergetic electrons moving towards 
the observer with the Lorentz factor 7, in general, has three power-law segments: a flat 
low-energy part, an intermediate-energy region which rises or slightly falls with a slope of 
less than unity (the exact value depending on 9), and the more steeply falling off part with 
the slope being between -2f3 and -2f3 + 1, again, depending on 9. 

The shape of the spectrum changes significantly with the angle 9 between the 
radiating particle's velocity and the axis of the current filamentation generated by the 
counterstreaming Weibel instability. As 9 — > 0, the low-frequency spectral break T\ 
approaches —00 and the maximum spectral slope (mid-range spectral index s i approach es 



the va lue of 1 (the trend of our results agreeing well with the 9 = case in 



Medvedev 



( 120061 )). As 9 increases, the spectral peak weakens as s\ decreases and T\ shifts towards the 



peak region. The disappearance of the spectral peak at some particular 9 appears to be a 
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result of both these spectral changes, and there is an extended transition region between the 
low-energy and high-energy power law trends. Consequently, we find that both the peaked 
and unpeaked spectra are well described by a three-region fit. Two-region fits are likely 
to miss out on the variation in the spectral slope below the peak at small 9; consequently 
the resulting low-energy spectral index will be extremely sensitive to where the peak falls 
relative to the lower bound of a measured spectral window. This will be true even if the 
low-energy spectral index is taken at a common energy, as in the "effective" low-energy 
spectral index a e ff commonly taken as the tangential slope of the logarithmic spectrum at 
25 keV. 

In comparing the radiation spectra for the full range of 9, we have found that the 
"drop point" , which we determined as the minimum of the numerical second derivative of 
the logarithmic data (i.e. the largest negative change in the spectral slope) serves as a 
good common reference point for both peaked and unpeaked spectra; in addition, the "drop 
point" in the radiation power spectrum evolves with the angle 9 much like the uF u spectral 
peak energy E p . 

In section [3j we developed in detail the relation between the radiation spectrum and 
the underlying Fourier spectrum of the particle's acceleration. In particular, we find that 
the radiation spectrum has much the same shape as the acceleration spectrum but that the 
apparent transition points in the two spectra do not simply coincide for most angles of 9. 
Furthermore, although the acceleration spectrum sees the re-emergence of a spectral peak 
for 9 — > 7r/2, the radiation spectrum does not. We have also demonstrated that a simple fit 
to the acceleration spectrum allows for the generation of a model radiation spectrum which 
approximates the realistic one with 10% accuracy. 

We have found that variations in the magnetic field spectral parameters influence the 
final radiation spectrum by controlling the width and peak-positions of the functions within 
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the integrand and the extent to which this directly modifies the effect of the offset, which is 
proportional to uj' . If we consider the general progression of the radiation spectra from being 
strongly peaked at small 6 to unpeaked at 6 near ir/2, the trends shown here indicate that 
the speed of the progression of the spectral shape between the two extremes is dependent 
on the relative strengths of the parameters in the magnetic field spectra transverse and 
parallel to the shock front. 

We have confirmed that the jitter radiation high-energy spectral index is determined 
primarily by the high-A; magnetic field spectral index /3, which otherwise has little influence 
on the spectrum. The low- A; magnetic field spectral index a is shown to have a significant 
influence on the low-energy and mid-range portions of the radiation spectrum when varied 
jointly in both magnetic field spectrum, although this influence is substantially reduced 
when only one a; L is varied. 

The parameters k± and k\\ represent the dimensionless correlation lengths of the 
magnetic field distribution in the direction along the Weibel current filaments (and the 
direction of shock propagation in the case of a GRB) and in the perpendicular plane 
(parallel to the shock plane for a GRB). We find that increasing k± and k\\ jointly shifts the 
entire spectrum to higher energies with relatively little effect on the spectral shape. Thus, 
as expected, the location of the spectral peak and break energies (and the corresponding 
peak energy E p of the vF v spectrum) are determined primarily by the correlation length of 
the magnetic field turbulence. The progression of the spectral shape between the head-on 
and edge-on cases is sensitive to the variation of the k parameter in one function relative 
to the other, such that for a particular viewing angle 9 either peaked or unpeaked spectra 
can be attained via modification of the k ratio K. In the extreme that k\\ is 2 orders of 
magnitude larger than k± we recover a peaked spectra for the angles as high as 6 ~ 80 
degrees. It is also notable that the spectral peak and transition points undergo relatively 
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little horizontal shift as K varies when 9 = 10 degrees, but shift quite dramatically (3-4 
orders of magnitude) during this variation for 9 = 80 degrees. 

We shall summarize the most notable properties of our jitter radiation results with 
particular significance to interpreting radiation spectra from astrophysical sources. First, 
the jitter radiation spectra are significantly harder than synchrotron spectra in the region 
just below the spectral peak. This may be a significant mechanism in astrophysical sources 
such as gamma-ray bursts where a substantial population is seen to violate the synchrotron 
limit. Second, both the maximum slope in the region below the peak and the extent of 
this sloped region (between the low-frequency spectral break and the peak) are strongly 
influenced by 9, with the largest slope and largest extent of the sloped region at small 9. 
The angle 9 np at which the peak disappears is determined primarily by the ratio between 
the magnetic field correlation lengths perpendicular to and along the filamentation axis. 
Third, the position of the spectral peak represents the characteristic correlation length 
of the magnetic field, which in case of the Weibel turbulence depends on the density as 
n 1 / 2 . This is at odds with the synchrotron radiation in which the spectral peak measures 
the magnetic field strength. This result is important for accurate interpretation of the 
observed spectra as well. Fourth, the high-energy part of the spectrum is represented by a 
power-law, even though the electrons are monoenergetic. Thus, no power-law distribution 
of Fermi-accelerated electrons is required to produce the observed power-law spectra in 
prompt GRBs. This has important implication for the interpretation of the observed data, 
provided the electrons are radiating in the jitter regime. Fifth, the angular dependence 
exhibited, combined with relativistic kinematics of a curved shock front can expla in certain 
puzzling features of the GRB prompt spectral variability (IMedvedev et al. 
sensitivity of the jitter spectra to the magnetic field anisotropy makes it a possible tool for 
diagnostics in sites of small-scale magnetic field turbulence and also a basis for analysis 
of astrophysical sources where the magnetic field orientations relative to the direction of 



200% . The 
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observation are changing over time. 



Although we often appeal to GRBs as sites where jitter radiation can likely be 
produced, we cannot exclude other astrophysical objects, e.g. jets in active galactic nuclei, 
early supernova shocks and other violent sources, provided that small-scale magnetic fields 
may be produced and maintained in them. At last but not least, one can use jitter radiation 



as diagnostic tool in 



2008 



Reynolds, et al. 



aser-p lasma interaction experiments, e.g., Hercules (IHuntington et al. 



20071 ). aimed at studies of Weibel turbulence and conditions in GRBs 



within the Laboratory Astrophysics and High-Energy-Density Physics programs. 
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Fig. 1(a). — A plot of the general form of the spectrum f(k) used in our paper for the 
magnetic field spectrum in the transverse and parallel directions. 
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Fig. 2(a). — The acceleration spectra ((w^j 2 ) for 6 ranging from 2 through 88 degrees. The 
spectra is numerically calculated for a step size of 0.1 in log(u/) and is normalized via division 
by the 3-dimensional integral over the magnetic field spectra. 



-34- 




0.5 



1.5 



2.5 



3 5 



Log fW'/Wn') 



Fig. 2(b). — View showing detail of peak region of acceleration spectra (jw^'l 2 ) for 9 shown 



in Figure 2(a) One can see the disappearance of the peak and flattening of the spectra for 
mid-range 9, followed by its reappearance at slightly lower log(u/) at 9 of about 76 degrees. 
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Fig. 3. — Plot of the spectral amplitudes taken at the low-u/ end of the calculated spectra, 
and at the approximate locations of the spectral peak for low-# and high-#. The dominance 
of one amplitude over the other illustrates the transition of the spectra from peaked to 
unpeaked as 6 progresses from to n/2. 
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Fig. 4. — The slope of the spectra in the log — log plot in Figure 2(a) Note that even 
for unpeaked spectra there is a flattening of the spectral slope and in some cases a local 
maximum (around log(u f / u f a ) = 2.5). 
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Fig. 5. — Plots of the product f(k)f(k + a) of two functions of the form of our magnetic 
field spectra (Equations ((Tj) and (J2J)), illustrating the effect of the offset a. As described in 
section [21 the acceleration spectrum is the integral of a product of such functions, with an 
offset controlled by the acceleration Fourier frequency uj' . 
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Fig. 6. — Here we show a model acceleration spectrum, as described in Equations (TTTjl - ( fT9l) . 
The spectrum is flat in Region I (log(cj') < ri), then becomes a power law u' Sl in Region II 
(ri < log(u/) < r 2 ), and a power law a/" 52 in Region III (log(u/) > r 2 ). 
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Fig. 7. — Here we show a comparison of our original calculated spectra (as presented in 



Figure 2(a)) for 0=10 (peaked) and 0=60 (unpeaked), and our model acceleration spectra as 
described in Equations [T7TfT9l In Region I, we extrapolate the initial calculated amplitude. 
In Region II, we fit the calculated spectra using the maximum slope in the region and fitting 
through the point of maximum slope. In Region III we have done a simple fit using the 
slope and position at point co>=4.6. Our choice of fit overestimates the peak of the calculated 
acceleration spectra. 
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Fig. 8(a). — The radiation spectrum for #=10 degrees as obtained analytically using equa- 
tions [20] -[22] and the fit values for the model acceleration spectra as shown in Figure [6] The 



result calculated (black) via full double numerical integration (as in Figure 9(a)) is shown 
here for comparison. The spectral transition point defined by equation [57] and the peak 
defined by equatiorj25l are marked with arrows. We note that the boundaries (dotted lines) 
between Regions I, II, and III as defined by our model acceleration spectra do not clearly 
correspond to transition points in the resulting radiation spectra. We have extrapolated 
slopes for Regions I, II, and III to demonstrate the variation. The model again slightly 
overestimates the peak energy, but matches the spectral shape well and agrees with the 
integrated spectra within about 10%. 
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Fig. 8(b). — The radiation spectrum for #=60 degrees as obtained analytically from our 
simplified model acceleration spectrum (as shown in Figure E]) versus the full double numer- 
ical integration result (black). As in the previous figure, we have extrapolated the slope in 
Regions I, II, and III. We again obtain a good agreement with the overall spectral shape, 
but now with an unpeaked form appropriate to this range of theta. The spectral transition 
point defined by equation [27| is marked with an arrow. 
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Fig. 9(a). — The angle-averaged radiation power spectra (F v ) of jitter radiation, numerically 
calculated for a step size of 0.2 in log(u;) for every two degrees in 9. The dotted lines are 
numerical calculations for values close to 9=0, illustrating the behavior of the spectrum in 
this limit. The spectra are arbitrarily normalized so that the first calculated value for the 
9 = 2 spectrum is 0. The radiation spectra are flat for low-cu, then slope upwards to a 
notable peak for low 9. As 9 progresses to 7r/2, the low-u/ amplitude increases and the size 
of the sloped region decreases until the spectra becomes relatively flat until its sharp decline 
at high-a;. 
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Fig. 9(b). — A more detailed view of the radiation spectra in the region of the spectral peak. 
The spectral values are calculated every 0.05 in log(o;) for finer resolution of the detail in 
this region. 
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Fig. 10. — The vF v spectra calculated from our jitter radiation power spectral results. The 
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Fig. 11. — A plot of the slope of the radiation spectra taken at several values of \og(u / 'co>o7 2 ) 
over the range of 9. 
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Fig. 12. — The angular dependence of the low-frequency amplitude A of our calculated jitter 
radiation spectra and the corresponding variation of the low-frequency amplitude A of our 
calculated acceleration spectra. The low-energy amplitude in each case is taken to be the 
first calculated value of the spectrum; for our choice of parameters and calculation window, 
this initial value is well below the first spectral break for all 9 > 2°. In both cases, we have 
normalized our spectra such that the low-energy amplitude at 9 = 2° is 0. 
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Fig. 13. — The angular dependence of the mid-range spectral indices Si of our calculated 
jitter radiation spectra and Si of our calculated acceleration spectra. For peaked spectra 
the mid-range spectral index is the maximum slope below the spectral peak; for unpeaked 
spectra the mid-range spectral index is the average slope between the point at which the 
spectrum falls below A — 0.01 and the point at which the numerical second derivative reaches 
its minimum value (the "drop point"). 
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Fig. 14. — The angular dependence of the high-frequency spectral indices S2 for our calculated 
jitter radiation spectra and 52 of the corresponding acceleration spectra. The high-frequency 
spectral index is calculated as the slope between the drop point (the position of the largest 
negative change in slope) and the higher-frequency position at which the numerical second 
derivative is closest to 0. 
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Fig. 15. — The angular dependence of the first spectral break positions \og e (u>/u 'y 2 ) = t\ 
and \og e (u)' / u' ) = T x in our calculated jitter radiation spectra and acceleration spectra, 
respectively. In each case these transition points are found as the intersection between the 
low-frequency fit line of slope and the mid-range fit line of slope s\ (radiation) or Si 
(acceleration), as found by the fit described in detail in section HI 
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Fig. 16. — The angular dependence of the second spectral break r 2 of the jitter radiation 
spectra and of the acceleration spectra. The transition point is the calculated intersection 
between the mid-range fit line of slope S\ (radiation) or S\ (acceleration) and the high- 
frequency fit line of slope — S2 (radiation) or — s 2 ) (acceleration). 
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Fig. 17. — A comparison of the angular dependence of the positions of the radiation and ac- 
celeration spectral peaks and the spectral drop points, where the numerical second derivatives 
of our calculated spectra reach a minimum. We have also plotted the angular dependence 
of the peak in the vF v spectrum. The drop point in our F v radiation spectrum nicely tracks 
the behavior of th e peak in the vF v spectrum, which is peak energy E p in the Band function 
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19931 ) commonly used to fit GRB spectra. 
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Fig. 18. — The angular dependence of the strength of the spectral peak, i.e. the height of the 
peak above the initial low-frequency spectral amplitude A. We see that the peak disappears 
in the radiation spectrum at 9 roughly 8° less than in the acceleration spectrum, and does 
not reappear at 6 close to 90°. 
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Fig. 19(a). — This figure shows the variation of the low- frequency amplitude of the radiation 
spectrum (vertical axis) with changes in the magnetic field parameters 014, which are varied 
progressively in the horizontal dimension as described below, for each of the fixed viewing 
angles 9 = 10°, 45°, and 80° (first, second, and third clusters respectively). The solid line 
with square data points indicates the behavior when a is jointly varied (over a range from 1 
to 10, as indicated) in both the transverse and parallel magnetic field equation: a = a± = an. 
The dotted lines show the effect of individually varying a± (triangular data points) and a\\ 
(diamond data points), which chracterize the magnetic field transverse to and along the 
instability filamentation axis, respectively. The parameters a± and an are varied from a 
starting value of 1.5 to an ending value of 2.5, in increments of 0.1. Our original value of 
a = aj_ = an = 2.0 is the central data point where the lines intersect. 
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Fig. 19(b). — This figure shows the variation of the low-frequency amplitude of the radiation 
spectrum (vertical axis) with changes in the magnetic field parameters which are varied 
progressively in the horizontal dimension for the fixed viewing angles 9 = 10°, 45°, and 
80° (first, second, and third clusters respectively). The solid line with square data points 
indicates the behavior when ft is jointly varied from 1.1 to 1.9, in increments of 0.1, in both 
the transverse and parallel magnetic field equation: ft = ft± = ft\\ . The dotted lines show the 
effect of varying ft± (triangular data points) and ft\\ (diamond data points) individually, also 
from 1.1 to 1.9, in increments of 0.1. Our original value of ft — ft± — ft» = 1.5 is the central 
data point where the lines intersect. 
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Fig. 19(c). — This figure shows the variation of the low- frequency amplitude of the radiation 
spectrum with changes in the magnetic field parameters k, which are varied progressively in 
the horizontal dimension for the fixed viewing angles 9 = 10°, 45°, and 80° (first, second, 
and third clusters respectively). We vary k jointly (k = k± = Kn) by powers of 10, from 1 
to 100, as indicated. We also vary k± and reii relative to one another by changing the ratio 
K = k±/k\\ through a range of values as indicated. 
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Fig. 20. — This figure compares the influences of the magnetic field spectral parameter varia- 
tions on the low-frequency amplitude A of the radiation spectrum obtained for representative 
viewing angles 9 = 10°, 45°, and 80°. For each spectral parameter (indicated on the bot- 
tom axis) the graph indicates the range between the maximum and minimum values of A 
obtained by our variations of that parameter. (The parameter variations are as indicated in 
the previous figures and described in detail in sectional) 
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Fig. 21(a). — This figure shows the variation of the mid-range spectral index (as determined 
by the fit described in section [5]) with changes in the magnetic field parameters aj. This is 
the maximum slope below the peak for peaked spectra, and the average slope in intermediary 
region for unpeaked spectra, and is shown for representative viewing angles 9 = 10°, 45°, and 
80°. The solid line with square data points indicates the behavior when a is jointly varied 
(over a range from 1 to 10, as indicated) in both the transverse and parallel magnetic field 
equation: a = a± = a\\. The dotted lines show the effect of varying a± (triangular data 
points) and ay (diamond data points) individually, from 1.5 to 2.5, in increments of 0.1. 
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Fig. 21(b). — This figure shows the variation of the mid-range spectral index with changes 
in the magnetic field parameters The solid line with square data points indicates the 
behavior when f3 is jointly varied from 1.1 to 1.9, in increments of 0.1, in both the transverse 
and parallel magnetic field equation: (3 = f3± = f3\\. The dotted lines show the effect of 
varying j3± (triangular data points) and f3\\ (diamond data points) individually, also from 1.1 
to 1.9, in increments of 0.1. 
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Fig. 21(c). — This figure shows the variation of the mid-range spectral index with changes 
in the magnetic field parameters k. We vary k jointly (k = kj_ = Kn) by powers of 10, from 
1 to 100, as indicated by the dotted line. We also vary k± and k\\ relative to one another by 
changing the ratio K = k±/k\\ through a range of values as indicated by the solid line. 
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Fig. 22. — This figure compares the influences of different magnetic field spectral parameters 
on the mid-range spectral index si of the radiation spectrum obtained for representative 
viewing angles 6 = 10°, 45°, and 80°. For each spectral parameter (indicated on the bottom 
axis) the graph indicates the range between the maximum and minimum values of s\ obtained 
by our variations of that parameter (The parameter variations are as indicated in the previous 
figure and described in detail in sectional) 
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Fig. 23(a). — This figure shows the variation of the high-frequency spectral index (as deter- 
mined by the fit described in section \5§ with changes in the magnetic field parameters CKj. 
The solid line with square data points indicates the behavior when a is jointly varied (over a 
range from 1 to 10, as indicated) in both the transverse and parallel magnetic field equation: 
a = a± = a\\. The dotted lines show the effect of varying a± (triangular data points) and 
an (diamond data points) individually, from 1.5 to 2.5, in increments of 0.1. 
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Fig. 23(b). — This figure shows the variation of the high-frequency spectral index with 
changes in the magnetic field parameters The solid line with square data points indicates 
the behavior when (3 is jointly varied from 1.1 to 1.9, in increments of 0.1, in both the 
transverse and parallel magnetic field equation: (3 = /3± = /3\\. The dotted lines show the 
effect of varying j3± (triangular data points) and j3\\ (diamond data points) individually, also 
from 1.1 to 1.9, in increments of 0.1. 
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Fig. 23(c). — This figure shows the variation of the high-frequency spectral index with 
changes in the magnetic field parameters k. We vary k jointly (k = = k\\) by powers of 
10, from 1 to 100, as indicated. We also vary and k\\ relative to one another by changing 
the ratio K = k±/k\\ through a range of values as indicated. 
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Fig. 24. — This figure compares the influences of different magnetic field spectral parameters 
on the high-frequency spectral index — S2 of the radiation spectrum. For each spectral pa- 
rameter (indicated on the bottom axis) the graph indicates the range between the maximum 
and minimum values of — S2 obtained by our variations of that parameter (The parameter 
variations are as indicated in the previous figures and described in detail in section O) 
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Fig. 25(a). — This figure shows the variation of the first spectral transition point (as described 
in section [S]) with changes in the magnetic field parameters «j, for 9 = 10°, 45°, and 80°. 
The solid line with square data points indicates the behavior when a is jointly varied (over a 
range from 1 to 10, as indicated) in both the transverse and parallel magnetic field equation: 
a = a± = a\\. The dotted lines show the effect of varying a± (triangular data points) and 
an (diamond data points) individually, from 1.5 to 2.5, in increments of 0.1. 
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Fig. 25(b). — This figure shows the variation of the first spectral transition point with changes 
in the magnetic field parameters Pi, for 9 = 10°, 45°, and 80°. The solid line with square 
data points indicates the behavior when P is jointly varied from 1.1 to 1.9, in increments of 
0.1, in both the transverse and parallel magnetic field equation: P = P± = P\\- The dotted 
lines show the effect of varying P± (triangular data points) and P\\ (diamond data points) 
individually, also from 1.1 to 1.9, in increments of 0.1. 



-67- 




Fig. 25(c). — This figure shows the variation of the first spectral transition point with 
changes in the magnetic field parameters for 9 = 10°, 45°, and 80°. We vary k jointly 
(/t = k,_i_ = K\\) by powers of 10, from 1 to 100, as indicated by the dotted lines. We also 
vary k± and k\\ relative to one another by changing the ratio K = k±/k« through a range of 
values as indicated by the solid lines. 
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Fig. 26. — This figure compares the influences of different magnetic field spectral parameters 
on the first spectral break log e (c<j/c<j 7 2 ) = T\ of the radiation spectrum obtained for repre- 
sentative viewing angles 9 = 10°, 45°, and 80°. For each spectral parameter (indicated on the 
bottom axis) the graph indicates the range between the maximum and minimum values of 
7~i obtained by our variations of that parameter. (The parameter variations are as indicated 
in the previous figure and described in detail in sectional) 
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Fig. 27(a). — This figure shows the variation of the second spectral break (as determined 
by the fit described in section ED with changes in the magnetic field parameters c^. The 
solid line with square data points indicates the behavior when a is jointly varied (over a 
range from 1 to 10, as indicated) in both the transverse and parallel magnetic field equation: 
a = a± = a\\. The dotted lines show the effect of varying a± (triangular data points) and 
an (diamond data points) individually, from 1.5 to 2.5, in increments of 0.1. 



-70- 




Fig. 27(b). — This figure shows the variation of the second spectral break with changes in the 
magnetic field parameters The solid line with square data points indicates the behavior 
when (3 is jointly varied from 1.1 to 1.9, in increments of 0.1, in both the transverse and 
parallel magnetic field equation: (3 = (3± = /3m. The dotted lines show the effect of varying 
(3± (triangular data points) and (3\\ (diamond data points) individually, also from 1.1 to 1.9, 
in increments of 0.1. 
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Fig. 27(c). — This figure shows the variation of the second spectral break with changes in 
the magnetic field parameters k. We vary k jointly (« = k± = k\\) by powers of 10, from 1 
to 100, as indicated by the dotted line. We also vary k± and K\\ relative to one another by 
changing the ratio K = k±/k\\ through a range of values as indicated by the solid line. 
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Fig. 28(c). — This figure compares the influences of different magnetic field spectral param- 
eters on the second spectral break log e (u)/uj 'y 2 ) = r 2 of the radiation spectrum. For each 
spectral parameter (indicated on the bottom axis) the graph indicates the range between 
the maximum and minimum values of r 2 obtained by our variations of that parameter (The 
parameter variations are as indicated in the previous figure and described in detail in section 

0) 
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Fig. 29(a). — This figure shows the variation of peak strength, namely the peak height 
relative to the initial amplitude, for changes in the magnetic field parameters «j. We present 
results only for 9 = 10° because at 9 = 45° and 9 = 80° the spectrum is unpeaked. The 
solid line with square data points indicates the behavior when a is jointly varied (over a 
range from 1 to 10, as indicated) in both the transverse and parallel magnetic field equation: 
a = a± = a\\. The dotted lines show the effect of varying a± (triangular data points) and 
ay (diamond data points) individually, from 1.5 to 2.5, in increments of 0.1. 
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Fig. 29(b). — This figure shows the variation of peak strength (the peak height relative 
to the initial amplitude) with changes in the magnetic field parameters The solid line 
with square data points indicates the behavior when j3 is jointly varied from 1.1 to 1.9, in 
increments of 0.1, in both the transverse and parallel magnetic field equation: f3 — /3± — /3\\. 
The dotted lines show the effect of varying j3±_ (triangular data points) and f3» (diamond 
data points) individually, also from 1.1 to 1.9, in increments of 0.1. 
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Fig. 29(c). — This figure shows the variation of peak strength (the peak height relative to 
the initial amplitude) with changes in the magnetic field parameters k. We vary k jointly 
(k = Kj_ = k») by powers of 10, from 1 to 100, as indicated. We also vary n± and k\\ relative 
to one another by changing the ratio K = Kj_/k» through a range of values as indicated. 
Note that we obtain peaked forms of the spectra for larger 6 when we have small values of 
the ratio K. 



